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We study the magnetoplasmon collective mode excitations of integer quantum Hall systems in 
a parabolically confined quantum well nanostructure in the presence of a tilted magnetic field by 
using the time-dependent Hartree-Fock approximation. For even integer filling, we find that the 
dispersion of a spin density mode has a magneto-roton minimum at finite wavevectors, at a few 
times lO'' cm~^ for parallel fields of order 1-10 Tesla, only in the direction perpendicular to the in- 
plane magnetic field, while the mode energy increases monotonously with wavevector parallel to the 
in-plane magnetic field. When the in-plane magnetic field is strong enough (well above 10 Tesla) ,we 
speculate that this roton minimum may reach zero energy, suggesting a possible second order phase 
transition to a state with broken translational and spin symmetries. We discuss the possibility for 
observing such parallel field-induced quantum phase transitions. We also derive an expression for 
the dielectric function within the time-dependent Hartree-Fock approximation and include screening 
effects in our magnetoplasmon calculation. We discuss several exotic symmetry-broken phases that 
may be stable in finite parallel fields, and propose that the transport anisotropy, observed recently 
in parallel field experiments, may be due to the formation of a skyrmion stripe phase predicted in 
our theory. Our predicted anisotropic finite wavevector suppression, perhaps even a mode-softening 
leading to the quantum phase transition to the anisotropic phase, in the collective spin excitation 
mode of the wide well system in the direction transverse to the applied parallel magnetic field should 
be directly experimentally observable via the inelastic light scattering spectroscopy. 

PACS numbers: 73.43.-f, 73.43.Nq, 73.43. Lp, 73.43. Cd 
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I. INTRODUCTION 



Observations of integral quantum Hall effect (IQHE, 1980) and fractional quantum Hall effect (FQHE, 1982) 
are important landmarks of condensed matter physics in recent decades [|l|. In quantum Hall systems, electrons are 
"frozen" (in their orbital motion) in discrete Landau levels by the external magnetic field, and have gapped excitations 
at integer or fractional filling factors. There is considerable richness of the phase diagram when additional (i.e. in 
addition to the orbital motion) degrees of freedom associated with spin, layer, or subband index are introduced [Q-^. 
These multicomponent quantum Hall systems have been extensively studied both theoretically and experimentally in 
recent years. In general, since the spin (Zeeman) energy is much smaller than the cyclotron energy due to the small 
effective g-factor and the small effective mass of electrons in GaAs based QH systems, the spin degree of freedom 
is not important energetically compared to the orbital motion. But spin can be crucial when a second quantum 
Hall system is coupled coherently (for example, in a double quantum well [DQW] system) or an additional magnetic 
field is applied in the direction parallel to the two-dimensional (2D) semiconductor quantum well plane. In the 
first situation, the finite barrier energy between the two wells opens a gap (As^s) between a symmetric and an 
antisymmetric subbands, which can be tuned by electron tunneling, layer separation, and/or bias voltage When 
^SAS is close to the Zeeman splitting energy, interesting physics has been predicted theoretically Q and observed 
experimentally H,^. On the other hand, physics of the second situation, where a tilted magnetic field is apphed 
to a wide width well (WWW) system to couple subbands of a wide well with spin-split Landau levels, has not yet 
been extensively explored. One reason for this is that the strength of the applied tilted magnetic field has to be 
very large (> 25 Tesla) in order to sufficiently enhance the Zeeman energy to be comparable to the Landau level 
separation in GaAs. Such strong and uniform magnetic fields has only been available very recently ||^. From a 
theoretical point of view, studying QH effects in a WWW with tilted magnetic field is difficult because the in-plane 
magnetic field hybridizes the 2D electron subbands arising from the confinement potential in the growth direction 
(i.e. perpendicular to the 2D plane) with the orbital Landau levels so that the electron wavefunction of a WWW is 
a complicated combination of electric ( "subbands" ) and magnetic ( "Landau levels" ) quantization even at the single 
particle level. It is sometimes simplistically believed that, if parameters are chosen properly in an isospin language, 
then a WWW system in a tilted field (at least for) the closest two Landau levels near the degeneracy point could be 
approximately mapped onto a DQW system. We emphasize that this mapping is not exact and misses subtle and 
interesting physics associated with a WWW in a tilted field. For example, experimentally a WWW in a tilted field is 
found to display both three-dimensional (3D) and two-dimensional properties ||]. In some situations a WWW system 
could behave very much like a DQW system albeit with strong tunneling) More strikingly, the recent observation 
of anisotropic resistance at even filling factors in a WWW system with an in-plane field M shows a possible stripe 
phase formation induced by electron-electron interaction near a degeneracy or a level crossing point. 

Inspired by the observed anisotropic transport properties at integer filling factors Q, we investigate in this paper 
the collective mode excitations of integer quantum Hall systems in a wide quantum well with a tilted magnetic 
field (i.e. in the presence of an in-plane magnetic field) by using the time-dependent Hartree-Fock approximation 
(TDHFA). We extend the work of Kallin and Halperin ||l^ for a strictly 2D system, i.e. a zero width well (ZWW) 
to a WWW system and derive a full analytical expression for the mode dispersion energy. To keep our theory 
analytically tractable we choose our quantum well confinement potential to be parabolic ("parabolic well"). Our 
choice of parabolic confinement is dictated by the fact that the corresponding single-particle problem (i.e. an electron 
moving in a one-dimensional parabolic potential along the z-direction in the presence of an arbitrary magnetic field) 
can be exactly analytically solved enabling essentially a complete analytic solution of the many-body TDHF solution of 
the collective mode spectra (essentially on the same footing as the 2D Kallin-Halperin work in Ref. ||l^]) in the WWW 
system in the presence of a tilted magnetic field (i.e. both the in-plane field and the perpendicular field producing the 
Landau quantization). The work presented in this paper is therefore a direct (and highly non-trivial) generalization 
of the strictly 2D Kallin-Halperin work ||l^ on the magnetoplasmons of a 2D electron gas (in the presence of only 
a perpendicular magnetic field) to a parabolic WWW system in the presence of a tilted magnetic field. We study 
both charge and spin mode collective excitations in systems of different electron densities, magnetic field strengths, 
and well widths. At even integer factors, we find that the dispersion of spin density mode has a magneto-roton 
minimum only at a finite wavevector in the direction perpendicular to the in-plane magnetic field, while it increases 
monotonously with respect to the wavevector parallel to the in-plane magnetic field. When the in-plane magnetic 
field is sufficiently strong, this roton minimum may reach zero energy before the ground state becomes polarized, 
suggesting a possible second order phase transition to a state with broken translational and spin symmetries. The 
possibility of this quantum phase transition (to an anisotropic symmetry-broken state) in the presence of a tilted field 
is one main new result of our work. We also derive the full formula for the dielectric function of the system within 
TDHFA by including the ladder diagrams consistently, so that it can be applied to other systems even when only few 
Landau levels are occupied. We include such screening in our collective mode calculation and discuss its effect to the 
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magneto-roton minimum. 

Before jumping into the details of the collective mode calculation, it is instructive to discuss in the appropriate 
context some earlier work in parabolic wells and in the ground state instability (i.e. the softening of collective modes) 
of similar systems. Among the models of finite width wells, parabolic wells are considered special, because the electron 
gas, in screening the parabolic conduction band edge potential, forms a constant density slab, being a good approx- 
imation to a 3D jellium where electrons move in a constant positive background charge density Furthermore, 
the parabolic confinement potential can be exactly diagonalized in a center-of-mass coordinate and therefore gives 
a non-spin-flip optical absorption energy exactly the same as its noninteracting result in the long wavelength limit 
(the so called generalized Kohn's theorem) [ p^p^ . As mentioned above,we use a parabolic confinement potential, 
because it allows us to find simple noninteracting eigenstates in the presence of a tilted magnetic field, which then 
provides a good starting point to consider many-body effects. The effects of imperfect parabolic confinement potential 
on the collective excitations have earlier been studied either with only a perpendicular magnetic field jl^ or with 
only an in-plane magnetic field [ p^ . Only rather small quantitative corrections were found (for example, small shift 
of resonance energy, and slight broadening of the absorption peak) for realistic wells (which necessarily deviate from 
ideal parabolic confinement considered in our work). We believe, therefore, that our theoretical results should a ppl y 
with quantitative accuracy to realistic parabolic quantum wells, and qualitatively to rectangular quantum wells pq ]. 

It is generally believed that in both three and two dimensions, when an infinitely strong magnetic field is applied, 
electrons undergoes a phase transition to a Wigner crystal state with broken translational symmetry at low tem- 
peratures. In the intermediate magnetic field region, Celli and Mermin [0 proposed a long time ago a possible 
exchange induced spin-density-wave (SDW) instability in a three-dimensional electron system. More recently, GaAs 
based semiconductor wide parabolic wells have been proposed as good candidates for observing such SDW instabil- 
ities since wide parabolic wells are essentially ideal 3D electron systems P Jl^ , p^ . Brey and Halperin proposed 
that the SDW instability and the transport anisotropy should be observed in a wide parabolic semiconductor quan- 
tum well system when an intermediate in-plane magnetic field is applied. Similarly, correlation-driven intersubband 
SDW instability has been predicted by Das Sarma and Tamborenea in DQW systems at low carrier densities [ pO[ . 
Intersubband-induced charge-density-wave (CDW) instability in a wide parabolic well with a perpendicular magnetic 
field was also investigated ||2l[]. To the best of our knowledge, however, these theoretically proposed (translational 
symmetry breaking) instabilities have not yet been observed experimentally. The only two experimentally observed 
candidates for charge (or spin) density wave instability in a quantum Hall system are the stripe phases (and the 
associate liquid crystal phases [^) in high half-odd-integer quantum Hall systems (i^ = 9/2, 11/2, etc) with or 
without in-plane magnetic field, and the stripe phases observed in an integer quantum Hall system in a wide well sub- 
ject to a strong tilted magnetic field 0. Although the ground state of the former system has been extensively studied 
p2| and is generally believed to be a "unidirectional coherent charge density wave" [^,^, the transport anisotropy 
in the wide well with a tilted magnetic field is not yet understood and not much theoretical work has appeared 
on this problem except for our recent short communication j26j. Our recent work based on Hartree-Fock (HF) 
calculation in a DQW system shows that spin-charge-texture (skyrmion) stripe could be the possible ground state for 
a WWW system, providing a possible explanation for the observed transport anisotropy in Ref. In this paper, 
for the first time we show the complete analytical and numerical work in calculating the collective magnetoplasmon 
mode dispersion within TDHFA and the observed mode softening confirms the existence of a novel phase proposed 
in Ref. g]. 

This paper is organized as follows. In Sec. || we obtain the single-electron eigenstat es in a parabolic confinement 



potential with a tilt ed m agnetic fi eld. We first discuss the noninteracting result in Sec. II A and then the interacting 



(HF) result in Sec. II B . In Sec. [I C we show that at even filling factors the system undergoes a first order phase 
transition from an unpolarized ground state for in-plane magnetic field, _B|| < B*^, where B*^ is a critical in-plane field 
strength, to a polarized ground state for _B|| > B*^. Based on the unpolarized integral quantum Hall ground state, the 



full theory with numerical results for the magnetoplasmon dispersion (within TDHFA) are given in Sec. HI. In Sec 



[V we derive the TDHF dynamical dielectric function for an integer quantum Hall system in a parabolic well with 
tilted magnetic field and use the result to study the magnetoplasmon dispersion in screened TDHFA. Implications of 
our results are discussed in Sec. ^ and finally we summarize our work in Sec. VI. 



II. SINGLE ELECTRON EIGENSTATES AND GROUND STATE ENERGY 

A. Non-interacting System 

We consider a parabolic confinement potential in z direction, Up{z) — ^ m*ojQz'^, where m* is the electron effective 
mass and luq is the confinement energy. A coordinate system is chosen such that the perpendicular magnetic field B± 
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is in z direction and the parallel magnetic field -By in x direction, with the 2D electron system being in the x — y plane. 

When the vector potential is chosen in a Landau gauge, A — (^,B\^x — i3||Z, 0), the noninteracting single electron 
Hamiltonian can be written as (we set ?i = 1 throughout this paper) 



Ho 
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eA 



2m* 
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+ l-m*ujQz'^ - g^iBBtotSz 
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(1) 



2m* 2m* 

where hb is the Bohr magneton, and g ^ 0.44 for GaAs. Sz is the z-component of the spin operator along the total 
magnetic field, whose magnitude is Btot — \J B\ + -B^- Py is a good quantum number in this gauge and can be 
replaced by a constant k (the guiding center coordinate). The remaining terms can be expressed by a 2 x 2 matrix 
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diagonalized by a canonical transformation, — U{0) ■ [x,z]^ and [px.Pz]^ — U{0) ■ [px.PzV ■> with 
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(2) 
can be 

(3) 



and tan(2^) = — 2a;_LW||/(w^ — lo\). The new Hamiltonian describes two decoupled one-dimensional (ID) simple 
harmonic oscillators in new coordinates, x and z: 



where 



Hq — 



'^1.2 



1 



2m' 



(pI+pI) 



m LOi _2 m LU2 _2 
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Using (n, fc, s) as eigenstate quantum numbers, where n — (ni, 712) is the orbital Landau level index and s ~ ±1/2 is 
the eigenvalues of Sz, one obtains the noninteracting eigenenergies, i?2 ^, and eigenfunctions, 0^ ^ s(^: 



(6) 



and 



aiky 



Ly S 



V'l^^^ (cos e{x + ilk) - sin Oz) ■ tp^^] (sin 9{x + l^k) + cos Oz), 



(7) 



= ^HAx + llk,z) 

where Ly is the system length in y direction and the function $^3(2; + l^k^z) has x and z components only. Iq = 
y^l/m*uj± = y^c/eB± is the conventional cyclotron radius. We keep the spin index in ^ because these notations 
will later be generalized to an interacting system, where explicit spin dependence may become crucial. In Eq. (|^), 
the function Tpn'^ (x) is defined to be 



with li = ^l/m*uji for i = 1, 2, and Iq = •\/l/m*w_L = a/ c/eB± is the conventional cyclotron radius. Hn{x) is Hermite 
polynomial. It is instructive to consider the asymptotic form of the eigenstate energies, Eq. and wavefunctions, 

Eq. (|^), in the following four extreme limits: (i) Taking an infinite well width limit, loq — > 0, wi — > Y/w^"^i-a^ and 
W2 ^ from Eq. Eq. (Q) then shows that the free moving direction is restored along the z direction, which is 
perpendicular to the total magnetic field, Btot, showing a 3D property, (ii) Taking a zero width limit (tJo oo), we 

have 9 i^l'i., cjq ^ oo, uj2 — > and therefore ■0tiV(^) ~^ and ^'^}{x) -tpnJix), the usual orbital 

wavefunction of a ID simple harmonic oscillator. Therefore by changing the value of ljq, one can obtain a quasi-2D 
system, which has both pure 2D and 3D properties by taking different limits of the confinement potential strength, 
(iii) Similarly, for zero in-plane magnetic field limit — > 0), we have uji —^ Max{uj±,LUQ) and uj2 — > Min{uj±,LUQ), so 
that the orbital motions in x and z direction are totally decoupled. This is the usual (i.e. without an in-plane field) 
quantum Hall system in a parabolic well, whose collective mode dispersion has been studied in the literature [ pT| . (iv) 
Finally we can take the strong parallel (in-plane) magnetic field limit (i3|| —^ oo), which is of interest in this paper. 
In this limit, we have 9 — > 7r/2, uji — > t^n — )■ oo, and L02 cjo^^i/^n ^ 0, i.e. the in-plane magnetic field enhances the 
effective confinement of a wide well system (compared to (ii)) and therefore a WWW system with a strong parallel 
field becomes similar to a thin well (strictly 2D) system with small Landau level energy separation. We emphasize, 
however, that our results shown below apply for any finite strength of valid to the lowest order of the ratio of the 
interaction strength to the noninteracting level separation. We will consider the strong in-plane magnetic field limit 
only when studying the screening effect in Sec. |V[ 

Energy levels described by Eq. (^ are shown in Fig. |l| as a function of in-plane magnetic field for a choice of 
parameters similar to the experimental samples in 0: electron density Ue = 0.42 x 10^^ cm~^, m* = 0.07 uiq (mg is 
the bare electron mass) and wq = 7 meV. The confinement energy is such that the size of the first subband electron 
wavefunction in zero field is 260 A. The perpendicular magnetic field, Bj_, is chosen to be 2.97 Tesla for v = 2. 

Using the noninteracting single particle wavefunction in Eq. (^, the noninteracting single electron Green's function 
can be easily obtained: 

G°(ri,Ti;f2,T2) = ^^G^,<^(ri,ri;r2,T2) 
- E E E 'l^lJir2)^lkA^i) [9iT2 - T^)9i^, - £;2 J - 9in - r2)0(£;2 , - ^)] , (9) 

(J k ^ 

where is the imaginary time, and /i is chemical potential at zero temperature. The Heaviside theta function 9{x) = 1 
for a: > and is zero otherwise. 



B. Interacting System in Hartree-Fock Approximation 

When electron-electron interaction is considered, we use self-consistent Hartree-Fock approximation (SCHFA) to 
calculate the single electron wavefunction self-consistently by including the Hartree and Fock potentials in the single 
particle Hamiltonian. This approximation is the standard leading-order many-body (self-consistent) expansion in the 
(unscreened) Coulomb interaction whose one- loop Feynman diagram representation is shown in Fig. ||(a). In 
SCHFA, the wave equation for the quantum Hall system is p8|: 



Ha + J df"V{r~- f") ^ l^m,p,s4'rn,p,J ir')4>rn,p,s{r') 4>fi,k,a{r) 
rh.p,s 

j dr'V{r- r')(l).n^k,a{r')^i^,%,p,a(l)m,p,aHr')<l>rn,p,cr{'r), 



(10) 



where Vm.p^s is the filhng factor at the specific quantum number, and it satisfies 



m,p,a 



(11) 



where is the total electron number. Tio is the same as in Eq. (|l|) by taking py = k. Note that the positive 
charge donor density (which produces the electron gas and thus provides charge neutrality for the whole system) is 
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not explicitly included above because these donors are usually located far away from the well in the experiment. In 
general, this background doping effect can be effectively included by introducing a screening length, A, into the bare 
Coulomb interaction, V{q), by writing V{q) = (47reVeo)(|9 P + {2Tr/Xf)-^^^. We take A = 620 A in our numerical 
calculation below to be comparable to the experimental setting [Q. This regularization of Coulomb interaction has 
little quantitative or qualitative effects on the results shown in this paper. The details of the donor screening and 
the exact value of A do not in any way affect any of our qualitative conclusions. For the situation wc focus in this 
paper, electrons are assumed to be uniformly distributed in the 2D well plane (i.e. i'm,p,s is independent of guiding 
center coordinate, p), and therefore Eq. (|l0|) can be simplified further as shown in Appendix To solve the SCHF 
equation, we first use the noninteracting wavefunction to calculate the HF matrix elements and then diagonalize it 
to get new eigenstates, which are used to calculate the HF matrix element again iteratively until self-consistency is 
achieved. The new single electron Green's function in SCHFA is similar to the noninteracting one in Eq. (^ except 
that the wavefunctions and energies correspond to the Hartree-Fock theory: 

Gn,a(ri,r2;ujj = — —j = }^<Pfi,k,a' [r2j(Pn,k,a[rij • (12) 

g° (uj) -E?^ k 

where ^2 (j(^) is Fourier transform of the noninteracting propagator, Q\ ^(O' 

Ql,{u^) = [ dte^-'Gl^it) = . i. _ (13) 



and similarly 



GnA^) = -■ ?rffF > (14) 



where the Hartree-Fock self-energy, = Ef; „ — , is obtained from the self-consistent solution of the Hartree- 
Fock problem (see Appendix in particular Eq. ( A2)). 



C. Level Crossing and Total Energy 

For our purposes, the most important feature of the spectra shown in Fig. |l| is the existence of a single particle level 
crossing at = 19.8 Tesla (for the noninteracting system), where toi ^ 36 meV, and u!2 — lo^ ^ 1 meV. The origin 
of this crossing can be understood by considering the asymptotic form of the energy levels in Eq. (j|) for large parallel 
magnetic fields: UJ2 ujqlo±/uj\\ — s- 0, so that at a critical in-plane magnetic field value (Sp U2 becomes smaller 
than the Zeeman energy leading to the level crossing shown in Fig. 0. For a non-interacting electron gas this level 
crossing leads necessarily to a (rather trivial) first order phase transition at _B|| = B*^ with an abrupt change in spin 
polarization for systems at even filling factors. Interesting quantum phase transition that may take place around this 
level crossing is the main subject of this article. In particular, we wish to investigate whether quantum level repulsion 
converts this first order transition to a second order quantum phase transition around this degeneracy point. Our 
calculated mode dispersion can be directly compared to inelastic light scattering spectroscopy, when such experiments 
are eventually carried out in these WWW systems in tilted fields. One can easily use our analytical results to study 
the magnetoplasmon mode dispersion at odd integer filling factors or for weaker in-plane magnetic field values (where 
intersubband coupling needs to be included). 

By using the self-energy obtained from Eq. (|l^) in SCHFA, the total energy of an interacting quantum Hall system 
can be obtained for a given electron configuration (7V| ,Ni), where iVo- is the orbital level index of the highest filled 
level of spin a. Considering double counting of interaction energy, the total energy in HF approximation is 



El <£;o. 



(15) 



To obtain the ground state energy, Eq, one should compare the total energies of all possible electron configurations 
and determine which one gives the lowest energy. As indicated in Fig. ^, a first order (noninteracting) phase transition 
from an unpolarized ground state (i.e. N-^ = = (0, 1^/2)) to a polarized ground state {Ni — N-^ + (0, 2)) is expected 
to happen at a critical in-plane magnetic field, BT,. In the third column of Table | we show our numerical calculation 
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results of B*^ obtained from Eq. (|T^) in the first order HF approximation for even filling factors = 6, 8. When 
the total electron density is fixed, is larger as the filling factor, v, is lowered by increasing the perpendicular 

magnetic field. Therefore our HF results qualitatively agree with the experimental data presented in Ref. except 
for a lower estimate of the critical magnetic field, B^^, which may be due to the correlation effects not included in the 
HF approximation and/or the nonparabolicity of the realistic confinement potential of the quantum well sample used 
in Ref. §. 

III. MAGNETOPLASMON EXCITATIONS 

In this section we will develop the full theory of magnetoplasmon excitations of an integer quantum Hall system 
confined in a parabolic well and subject to a tilted magnetic field within TDHFA. For zero width (pure 2D) wells with 



a perpendicular magnetic field only, magnetoplasmon modes were investigated in |10|. We note that magnetoplasmon 
excitations in parabolic wells have been theoretically discussed previously in the literature |f^jl9| , pl] , p9| -pi| in different 
limited conditions. Our work goes beyond results presented in those papers and we derive the exact dispersion of 
collective modes in the lowest order of the ratio of Coulomb interaction to the noninteracting Landau level separation. 
In a WWW with tilted magnetic field, there is no translational symmetry along the growth direction (z), which is 
hybridized with the in-plane components (x — y) so that a many-body theory developed in momentum space seems 



not to be particularly useful. However, it is shown below in Sec. [II A that the in-plane momentum of an electron-hole 
dipole in such WWW with tilted magnetic field is still conserved, showing the existence of a well-defined electron-hole 
bound state (a magnetic exciton) |I0[| and the collective mode dispersion along the 2D plane can still be obtained 
analytically as we show below . The full many- b ody theory and the numerical results for collective mode dispersion 
are shown in Sees. IIIB - |HI E and in Sec. HI F respectively. 



A. Momentum Conservation of an Electron-hole Dipole Pair 

As pointed out in Ref. jf^] , a crucial fact that allows one to explicitly write analytical expressions of the energy 
dispersion of magnetoplasmon excitations is the existence of a good quantum number in the problem given by the 
well defined in-plane momentum of the electron- hole dipole pair (magnetic exciton). It is easy to show that their 
argument can be extended to the case of a WWW with an arbitrary confinement potential along the z direction 
even in the presence of a tilted magnetic field. This is not obvious since the tilted magnetic field typically hybridizes 
the in-plane motion with the subband dynamics perpendicular to the plane, destroying the apparent translational 
symmetry. Consider the Hamiltonian of a magnetic exciton or an electron-hole pair in a general quasi-2D system, 



2m* 



(pi - -A{r,)f + (p2 + -A{r2)f - V{r, - r^) + C/(zi) + U{z^), (16) 



where particle momenta pi, vector potential A, and particle coordinates are all three dimensional vectors. V{r) 
and U{z) are electron-electron Coulomb interaction and the quantum well confinement potential respectively. The 
Zeeman term is neglected here because it is irrelevant for this discussion. Following [ p^ a magnetic exciton momentum 
operator can be defined to be 

Qx=Vi+P2~ -{A{ri) - A{r2)) + -Btot x (ri - ra), (17) 
c c 

where Btot — B^z + B\\x. Using the Landau gauge for the vector potential one can easily verify that the in-plane 
components Qx.i^ — {Qx,xiQx,y) commute with the Hamiltonian. Existence of dipole excitations with well defined 
momenta (eigenvalues of Qx,l) immediately follows from this commutation. Similar to Ref. [l0|] , we can construct 
the zero-momentum magnetic exciton wavefunction in a parabolic well with a tilted magnetic field: 

vI/^';;"£(Ax,A2/,Z,Az) = j d77e-*'''^^/'o<i>_^^^^^(^ + A2;/2,Z + Az/2)$^i„,,j77-Aa;/2,Z-Az/2), (18) 

where a hole is in a state Ha — n and a particle is in a state ftfj — fia + rfi. For the exciton wavefunction of finite 
momentum, one just needs to replace Arj^ — (Ax, Ay) by Arj_ — l^q^ x z, and introduce a plane wave prefactor 
for the center of mass coordinate [ pO{ . Note that this wavefunction has additional dynamics along z direction: center 
of mass {Z) and relative (Az) coordinates of the electron-hole pair. Other than this additional z dynamics, the only 
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difference between the exciton wavefunctions for the ZWW 2D system '10| with a perpendicular magnetic field (see 
Eq. (Bl) or jl^) and for the WWW with a tilted magnetic field (our interest in this paper) is that the latter has 
one more Landau level quantum number associated with the subband dynamics induced by the confinement energy 
of the well. In Appendix ^ wc show that the magnctoplasmon energy in our theory can be expressed in terms of 
the magnetic exciton wavefunction given in Eq. ([l^). This provides a more comprehensive and physical picture for 
understanding the collective mode excitations discussed in this paper (see Appendix § in this context). 



B. Correlation Function 



In the linear response theory, collective mode energies are obtained by the poles of a density correlation function, 
Hx, where A = p, S± and Sz, for the singlet charge density mode and the three triplet spin density modes respectively. 
We define an operator, @x = 1, 2S±, 2Sz, respectively for the spin vertex operator of each corresponding correlation 
function. In this notation the most general form of these correlation functions in coordinate space is 

Ux{f,t;f',t') = -z^^[e,].„.je,]..,., (t [*t^(f,t)^'^^(r,t)#t,(f',t')^^,(f',t')]^^, (19) 
0-1,2 (t;_2 

where 5'J.(r, i)(5'(j(r, t)) are the electron field creation (annihilation) operators of space r and spin a at time t; T[- ■ ■] 
is time-order operator, and (• • ■)c is the expectation value of the interacting ground state. In a WWW system, 
there is no translational symmetry along z direction so that one has no correlation function in momentum space in z 
direction. The usual momentum space description for the vertex function and the related Dyson's equation then seems 
not feasible because the in-plane magnetic field mixes the z dynamics with in- plane dynamics jl^,^,^ . Actually the 
system is more like a 2D quantum dot [p^ in x — z plane confined by two independent parabolic potentials along x 
and z axes as shown in Eq. (^. The method we develop in this paper, however, enables one to obtain directly the 
appropriate Dyson's equations for the screened interaction and the vertex function without evaluating the correlation 
functions of Eq. (|l9|). The magnctoplasmon excitation dispersion and the dielectric function relevant for screening 
can be read out directly from our equations given in the next section. Note that the theory developed below is 
independent of the exact form of the single electron wavefunction and is completely general within the TDHFA. 

C. Screened Interaction and Vertex Function 

Before exploring the many-body theory for the collective mode, we first define the interaction matrix element, which 
will be used frequently later. Using the interacting single particle wavefunction, the unscreened matrix element of a 
bare Coulomb interaction, V{f), can be obtained 

^HaMrs '''''''' ^'^^'^^^'^^''■^ ^ J '^^^ J ^^^2 V"(fi -r2)0Si,fci,o-i^(?'l)'/'n2,fc2,O2^(n2)0fi3,fc3,O3(n2)0«^ 

where f2 is the well volume and we define an effective interaction, y5i5^ ^ - (q) = V(q)AZ'-Z^ (—q)AZ?''-? (q) where 
the form factor AZ^Z' [q) is obtained from the single particle wavefunctions 



dx [ dze-'''-'-'-''>^'^n^,,^{x - qyll/2,z)'^fi„a,{x + qylll2,z). (21) 



Momentum and spin conservations during the scattering process have been included in Eq. (|2£ 

Note that in the strong parallel magnetic field regime, (10 Tesla < B\\ < 25 Tesla), only a few Landau levels of the 
first subband (rii — 0) are occupied at zero temperature since lui ^ W2 (see Eq. (^) and Fig. |l|). Therefore we could 
omit the first orbital (i.e. subband) level index and neglect all intersubband transitions (i.e. excitations between levels 
of different ni) by assuming for simplicity that rii = throughout our analysis and numerical calculations shown 
below except where noted otherwise. In other words, the vector representation used in Eq. (^) for orbital Landau 
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level index, n = {ni,n2) = (0,77.2), is simplified to be n and so are all other orbital notations (like m = (0,7712) = ttt, 
and N„ = {0,Na) = Ncr, etc.) from now on in this paper. It is straightforward to extend all of our analytical and 
numerical results to include excitations of both orbital quantum numbers. All analytical results would retain the same 
form with additional level indices (i.e. other value of rii) showing up in the formula. Our numerical results will not be 
affected at all by this assumption in the strong in-plane magnetic field region of our interest where hybridization with 
higher ni levels is negligiblly small. We also will not show the spin index explicitly during the derivation except in the 
final results. We start from the screened Coulomb interaction, V(fi,r2iti — t2), caused by electron-hole polarization 
(see Fig. 1(b)): 

Viri,r2;ti-t2) = Viri-r2)S{h-t2) + J dfs J dfiVin ~ r3)W3:ti;ri,t2)Vifi ^ r2), (22) 

where n(fi, ti; f2, t2) = np(fi, ii; r2, i2) is the reducible charge polarizability (see Eq. (p^)). Multiplying by single 
particle wavefunctions and doing the space integration, Eq. (p2|) can be transformed to 



Vi,4:2,3(ii - ^2) = J dri J df2V{ri,r2;ti ~ t2)4'i\ri)(f)2H'f'2)(t>3ir2)4'4.iri) 

= ^l,4;2,3'5(il - ^2) + ^ Vi^4;l3a / dtsdteGaih - t5)Q (3{tQ - ti) 

0« ''' (f^5 ) '/'/3 (ns )7 (^5 , ^5 ; ^3 , ^6 ; f^4 , ^2 ) ^(?^4 - 7^2 ) '/>2 ''' (7^2 )03 (7=*2 ) 



'2,4,5,6 

= "l^l,4;2,3'5(il - ^2) + ^ ViA-f3af^f3,2Atl " h), (23) 

where we have introduced a conventional reducible vertex function, 7(^5, t5; 7^6j ^e! ^2), in coordinate space to express 
the reducible polarizability, 11; index a{(3) denotes all related quantum numbers of that level, {ma(0),Pa{i3), o'a(f3)), i.e. 

Mr) = 0m„,p„,<T„(r'), Qait) = Gm^^a^it), 14,4;q/3 = Vnin^',mfm'^''''°', and faf3;2,3 = ^m^m^-nXs""'^ '''^'"^ for simplicity 

(number indices represent external variables, while Greek indices represent dummy variables in a summation) and 
= J dri. To avoid confusion, we clarify our notations which are necessarily different from the standard many-body 
textbook terminology because of the highly complicated nature of our single-particle wavefunctions. First Vi.4;2,3 and 
Kiin4,n2n3 (*?) different functions according to their definition in Eq. (|20|); secondly the f function in Eq. ( p3| ) is 
not the same as the conventional definition of a vertex function due to our inclusion in F of additional two electron 
Green's functions and one interaction term (however, 7 is the same as the conventional reducible vertex function in 
coordinate space. See Fi^ I (c)). This is because, unlike a ZWW (pure 2D) in Ref. Q or a WWW without any 
in-plane magnetic field |2l| , the z-component of electron wavefunction of our system is not separable and therefore 
cannot be ignored. It is more convenient to work in the relevant conserved quantum number space rather than in the 
conventional momentum space. 

The leading order (of the ratio of the interaction strength to the noninteracting energy separation, W2) of the vertex 
function, Tai3-2.3{ti - h), is obtained by using 7(r*5, is; r*6, te; ^4, ^2) = <5(^5 - ^4)<5(r6 - ri)S{t5 - t2)S{te - 12) in Eq. 

^af3;2.3{tl - ^2) = Ga{tl - t2)Q[3{t2 — ti)Vai3-2,3, (24) 

which has the following Fourier transform in time: 

rQ/3;2,3(w) = Da(3{l^)Vai3-2,3, (25) 

where (after retrieving the spin index) 

Daf)[uj) = -. ^_ ^ (26) 

is nonzero only when the dipole pair, [a, (3), represents one hole in the filled level and one electron in the empty level 
at zero temperature. To avoid confusion, here we clarify the meaning of crQ(^), rn^^^^^^ and N^^^p.^ in Eq. ( pq ) again: 
Cq(/3) is the spin quantum number of state a{(3), rria^^i^^ — (0, 777,ct„(^) ) = ™o-„o) is the orbital Landau level index of 
state a{(}), and N^^^^^ = (0, N^^^^^) = N^^^^^ is the orbital Landau level index of the highest filled level of spin crQ(0) 
as first defined in Eq. ([T5|). We will see later that Da[i{uj) is the only dynamical part of the vacuum electron-hole 
bubble. Note that when spin is included, V^p^jx implies CTq, = ap and = ax automatically because of the manifestly 
spin-conserving non-spin-fiip nature of Coulomb interaction. As a consequence, Ta/3-^2,3{ti — ^2) in Eq. ( |25| ) (but not 
7(7^5, ts; 7^6, ^e; ^4, ^2) in Eq. (psh) then becomes identically zero when considering spin-flip excitations. 
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D. Dyson's Equations in TDHFA 



Including ladder and bubble diagrams as shown in Fig. ||(c), Dyson's equation for the full vertex function, 

rQ/3;2,3(il - h), is 



rQ/3;2,3(il — ^2) — ra/3;2,3(^l — ^2) + / dt^dt^Qaiti — t5)Qfj{tQ — ti)S{t5 — tg) 

dtrdhG{b, 7)G(8, 6)y(f5 - r&)l{r7,tr, rg, ig; ^4, ia) 



X / 0a^(r5)0/3(n5)F(f4-r2)02^(f2)03(r2) 
"'2,4,5,6 



7,8 



X / 0a^(r5)0/3(r5)V^(r4 - r2)02^(?^2)(?!)3(r2) 

J2,4,5,6 



dtjdtsG{6, 7)G(8, 6)^(7% - r6)j{r7, ^8, ^8; ^4, ^2 



(27) 



(28) 



(29) 



(30) 



Eq. (|29| ) and Eq. (|30| ) are respectively Dyson's equations for the vertex function and the interaction matrix element 
in the quantum number, a = {ma,Pa,o'a), space. 

In order to investigate the magnetoplasmon dispersion, one has to integrate out the continuous variable, fc, in 
Eq. (^9|) and Eq. (^) to get a matrix representation in the level index only. Taking into account the momentum 
conservation shown in Eq. (|20|), we define a new unscreened matrix element and a new bare vertex function given by 
(let q± = qxX + qyjj be the in-plane momentum) 



which can be further simplified to 

f a/3;2,3(^l ^ ^2) = rtt^;2,3(^l — ^2) + / dt^Qaiti ~ t5)Q/3{t5 ~ ^l) ^ [^Kt^^^^ + Va{i-vii\ f piy;2,3(^5 — h) , 

with the following Fourier transform in time 

f a/3;2,3('J^) = ra/3:2,3(^) + Dap{uj) ^ [—Va^i-yp + Vap-vfi] f /ji/;2,3 (w) . 

Similarly Fourier transform of Eq. ( p3| ) gives 

^^l,4;2,3(^) — ^1,4;2,3 + ^ yi,A;Pa^ ap-2,^{^) ■ 

al3 



TJ = pl(kl-k2)qJi Yki+qy/2,ki-qy/2;k2-qy/'2M+qy/'2 — _}^_\^ Y Co, n' \ 

'^nirn,n2n3\M±) — / , ^ nirii, 112113 ~ 2'Kp L ' '^"l"4,"2"3 Vj-L ; 'Jz 7 



(31) 



I'z 



A, 



mo.in0;n2ii3 



(r7, ,,j\—\r) pKkl~k2)qJl Yki+qy/2.ki-qy/'2\k2-qy/2M+qv/'2 _ TT 

) — [J-^inainfi\yJ )\ ^ mc,mi3-n2n3 — <^ inc:,mfi -,112113 {q± ) , \'J'^) 



where L^Ly — 2TrlQN^, and is the degeneracy of each Landau level. Note that in Eqs. (|T]) and ( |3^ ) only the in- 
plane component of momentum, q± , is shown explicitly. This follows from the fact that the in-plane exciton momentum 
is a good quantum number even in the presence of tilted magnetic field as discussed in Sec. [II A. Expressions for the 
screened matrix element, U, and the full vertex function. A, can be similarly obtained using the fc-summation over 
V and r as in Eqs. (|l|) and After some tedious analysis, we obtain the following pair of matrix equations in 

Landau level indices 



111114,112113 



A. 



714,712 713 (^Ql. ) ^ ^ Unin^^m^ma (^-L ')-^mamj^ ('^) -^771^771^ ;7i2 7i3 (^J- 7 ^) ? 

TllaTTl^ 

aini3;in„mij {q±)D„i^„i^ (w)Am^Tn„;n2n3 



where the new interaction function, W, is 



(33) 
(34) 

(35) 
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The ladder (exciton binding) energy, 17^™'^ ((fj^), and the random-phase- approximation (RPA) energy, U{q±), are 
respectively (after retrieving the spin index) 

q' q' 

^ g/ ^ 

Note that the non-spin-flipping interaction, (ctq, — cr^, = (^p) for Eq. (|3^) and {a a = C/SjCiy — (Jp) for Eq. (|3^), is 
already incorporated above. As mentioned in Sec. Ill B| , instead of calculating the irreducible polarizability directly. 



we derive the Dyson's equations of the interaction matrix element and a special vertex function in Eqs. ( |33| ) and 
(^), which can be used to obtain the collective mode energy and dielectric function. Above derivation and results 
are independent of the details of single particle wavefunctions or eigenenergies, and are valid for arbitrary quantum 
well confinement potential, provided the relevant form function A^inj ('f) of Eq. ( |2l|) is appropriately modified. In 
considering the spin degree of freedom, only non-spin-flip modes are included in Eq. (^3[) and Eq. (^^. The spin-flip 
modes do not include the bubble diagram due to the spin conservation implied by the interaction shown in Eq. (p3). 



and Wyn^rni3-m„m^{(i±) of Eq. ( |35D therefore becomes the same as the exciton binding energy, ^m^m^l^Dj foi' 

the same reason. For convenience we will not distinguish these two modes here until we get to the final results in the 
following section. 

E. Energy Dispersion of Magnetoplasmon Excitations: analytical expression 

Solving Eq. ( ^ ) one can obtain the vertex function, A„i^„i^;„2„3 ((fx, lj), and substitute it in Eq. ( |3^ ) to get the 
full formula for screened Coulomb interaction (using Arn^m^:n2n3{^j-> , q±) = Umf^m^;n2n3{q±) according to Eq. (|32|)): 

Unini,n2n3 ('^j 9±) 

= CAiin4,n2n3 (q±) ^ ^ ^ Unin4;7nf^ma{Q-l-)-^'mc^nifj{^) ^ ^ [^q/i(^/3i/ -^rriam^ (^) ^^m^m,y ;m/3ma (q± )] A^t^^^^ ^^^ns (q_L i 



where the dielectric function, e{uj,q±), is a matrix function, 



and the "dispersion matrix", Y, is 

The TDHF dynamical dielectric function appearing in Eq. (^9|) and Eq. ( pO[ ) includes infinite series of both RPA 
bubble diagrams and the excitonic ladder diagrams. Theoretically given a finite matrix size by including relevant 
Landau levels (i.e. by appropriately cutting off the infinite matrix equations give above), one can numerically calculate 
each element of the dielectric function and obtain the collective mode dispersions by solving the standard collective 
mode equation, (ie<{e(w, ql)} = 0. However, it is easy to see from Eq. (|39|) that solving lo from det{e{uj , q^yi = is 
the same as solving uj from det{Y{uj, (f_L)} = 0, or more conveniently, the same as solving the eigenvalue equation of 
y(0, qj) because Y{uj, q±) = Y{0, q±) + iw ■ I, where / is the identity matrix due to the special form of D,n„mf) (^) in 
Eq. (Eq). Therefore focusing on the collective mode dispersion in this section, we will discuss the dispersion matrix 
i^(0, 9J-) below in more detail instead of the dielectric function itself, which is studied in the next section. We note that 
this theoretical simplification of the equivalence between the simple static Y function and the dynamical dielectric 
function in obtaining the collective mode dispersion has not earlier been appreciated in the literature. 

According to Eq. (|2^), the only valid matrix element of Eq. ( |40[ ) should be for the pair, (mc, mp), of one electron in 
an empty level, ma{mfj), and one hole in a filled level, m/3{ma). To the lowest order of (e^/eo^o)/w2, only four levels, 
(iV,t), (iV,i), (7V-M,t), and (7V-M,i) are included (see Fig. |(a)), where = Z//2-1 is the level index of the highest 
filled level (note the first Landau index has been taken to be zero). After separating spin-flip and non-spin-flip modes. 
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one can obtain two 2 x 2 Y matrices for spin-flip (a) and non-spin-flip (p) excitations respectively (after retrieving 
the spin index): 



jjbind,l] N 
'-'AT AT wj-iTV+l Vii / 

N+lNM+lN^'l^) 
bind, il 



-u 



A pTi 



bind, II 
NN+1,NN+ 



l(<ZJ 



bind, I 
NN,N+1N+ 



u 



NN,N+1N+ 



u 



n 



NN+1,N+1N 



NN+1,N+1N 



as" 



^AW + 1 , AT -I- 1 Af ( ) 



TT 

AfAf-|-l,Ar+lN 



where AE"^"^ — 



-E° 



^niai is the HF single particle energy difference. Note that the off-diagonal 



term in ^^(91) is omitted in the paper by Kallin and Halperin |^ for a ZWW system, because they just considered 



the leading order 1x1 matrix representation of Yo-. Using ?7^"'|['^^^_|_]^^((7j^) 



bind, 11 



NN+l,NN+li'lA- 



we can obtain 



the dispersions of three spin collective modes and one charge collective mode accordingly by solving the determinantal 
equation for Y. (Note that for systems of even flUing factor, v = 2{N + 1), the spin index can be neglected in the 
self-energy, S^^, and interactions, U and Jj''"'-'^^ if the ground state is unpolarized and hence spin-symmetric): 





= UJ2~ 


, yHF 


^pHF 


jjbind 
'^NN,N^ 


-17VH 




wct-('?±) 


= Ll>2 ^ 


, yHF 


-^HF 


J jbind 


-17VH 






= UJ2 ^ 
= L02 ~ 


, yHF 
, yHF 


•r^HF 

— 

•r^HF 


J jbind 
~ '-^NN,NA 

J jbind 
~ '-^NN,N^ 


-17VH 
-17VH 


Ml-) 



Tjbind (ff •• 



/4 



J jbind 



liqi 



/4 



2Unn 



+ l,Af + lW 



(43) 

(44) 

(45) 
(46) 



Equations (|4^)-(4^) are our main analytical results in this paper, which are the formal generalizations of the 
corresponding results in Ref. |l^ to a WWW in the tilted field, except for the square-root term in Eqs. ( ^ ) and (|44|). 
This shows that the off-diagonal term of Ya{q±) is an exchange-interaction- induced level-repulsing effects for the two 
spin triplet modes, uja-± , and effectively increases the Zeeman energy. Note that the analytical derivation given above 
is not related to any specific form of the single particle wavefunction which enters only through the actual calculations 
of the various matrix elements in Eqs. (^)-(^). The only constraint on the wavefunctions is that they must be 
obtained in a conserving approximations. There are several approximations we can use to obtain the single particle 
wavefunctions and eigenenergies. Here we will compare two of them: one is the fully SCHF approximation as shown in 
Eq. ([l0| ) , and the other one is the first order Hartree-Fock approximation, where electron Hartree and Fock potential 
are calculated by noninteracting electron wavefunctions, which are not renormalized by a self-consistent equation (see 
Fig. 11(d)). It is shown later that such a first order Hartree-Fock approximation does capture the most important 
contribution of the SCHFA, but is computationally much easier than the numerical results from the SCHF equations. 
In fact, as mentioned before, the TDHFA in solving collective mode dispersion is exact only to the leading order in the 
interaction. Therefore in some sense the single particle wavefunctions and eigenenergies calculated in the full SCHFA 
are not guaranteed to give better collective mode dispersion energies than those calculated in the simple first order 
HF approximation, although the former may very well be better in calculating the single electron properties (such as 
the electron density profile or absorption spectra [^). In fact, we believe that in the spirit of our TDHFA calculation 
for the collective mode dispersion, it is actually better to use the first order HF wavefunctions and energies in the 
collective mode calculation in view of the excitations of the theory in the leading order Coulomb interaction. The use 
of such first order HF wavefunctions and energies in the TDHFA calculation of collective mode dispersions ensures 
that all quantities entering the theory are leading order in the Coulomb interaction . Therefore it is instructive to 
show the corresponding formula of the first order HF approximation in our theory here and we will compare the two 
sets of numerical results (SCHF and first order HF) in the next section. Defining the first order interaction matrix 
element similar to Eqs. (gO|) and (|2l] ) by using the noninteracting wavefunctions in Eq. (|^), we have 



dft 



dx / dz e 



-9y/2,l 



'^fij ,qy /2.aj ) 

Xx~qyll/2,z)^%^ 



{x + qyll/2,z), 



(47) 



whose analytical expression for a parabolic well could be obtained by using the generalized Laguerre polynomial 
discussed in Appendix^. As a consequence, one can also obtain the analytical expression corresponding to Eqs. (^^- 
( p6| ) in the first order Hartree-Fock approximation. For convenience, we first define two new dimensionless quantities, 
Qi{q) and Q2{q), as following: 



12 



Qiiq) 
Q2{q) 



COS^ 9{qylQY + {cosOqxlo — siT\9qzh)'^\\ 

sin^ 9{qylo)'^ + (sin Oq^lg + cos 0qzlo)^^2 
2X0 ■ 



(48) 
(49) 



The first order RPA (direct) energy then becomes (suppressing the spin index here): 
1 I ^ 



= 2t^i2i 1) ^ [-^1(9-^' 9^)] exp hQ2('f±, 9^] Q2{q±,qz) \Llf {Q2{q±,qJ)\^ , (50) 

where LJ^(x) is the generalized Laguerre polynomial and the first order exciton binding (exchange) energy is 

U^N^J:Nt^N+liq±) = ^ E "^^^l-l'y ~ qyq^)ll)V{q')A%{-q')A^^l,^^,iq') 



= ^ E - qyq'.)ll)V{q') exp [-Qi(<fl)] cxp [-Q2(gl)] (Q2(gl)) i^V+i (^2(9!)) (51) 

?' 

For the first order HF self-energy, it is more convenient and instructive to show the self-energy difference between 
levels N and + 1 individually for the direct or the Hartree term (S^) and the exchange or the Fock term (S^): 



^N+l 



N 



■ E E ^(9^)4"^ (9.) A^iN+liq^) - A%iq.) 



1=0 



N 



J2 V{q.) exp i-QM)] exp [-Q2iq'J] [L%^^ {Q2{q'.)) - L% (g2(g;))] E ^? (QM) , (52) 



9^ 



/=0 



N 



N 



.1=0 



1=0 



-1 1 ^ 

— E Viq) exp hQiiq)] exp [-Q2{q)] ^ E(Q2(g))^-';! 



It is easy to prove that 



iV! 



^l-^Q^[q)\L^+^-\Q,{q))\^ 



\L 



N-l 



{Q2{q))\ 



(53) 



NN+1,N+1N (0) 
{l),bind 



-'N+1 — NN,N+lN+li^) 

by using the following two identities for the generalized Laguerre polynomials: 

xL^+\x) = (n + m + l)L^{x) -{71 + 1)L™, i(x) 



(54) 



(55) 



Eq. ( |54| ) shows that in the long wavelength limit, the charge density collective mode has the same energy as its 
noninteracting result (as it must). 



Lo^'Hq±^o) = K 



N+l,l 



(56) 



which reflects the generalized Kohn's theorem |1J] . This shows that the time-dependent Hartree- Fock approximation 
we apply in this paper is a current-conserving approximation to the leading order single electron wavefunctions and 
eigenenergies. From the numerical calculation presented in the next section, such a generalized Kohn's theorem, 
^p{q± ^ 0) = W2, is true also for Eq. (^), where the electron wavefunction is calculated self-consistently through 
Eq. (p^). However, one should note that if one includes the larger matrix size in Eq. (|4^) to go beyond the lowest 
order in (e^/eo^o)/i^2 (see Fig. ^(b)), there is no such exact cancelation, since some more diagrams (higher order in 
the interaction) should be included in Fig. |^ in order to obtain the current-conserving theory for collective modes in 
higher order calculations. 
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F. Energy Dispersions of Magnetoplasmon Excitations: numerical results 



The two 2x2 matrices shown in Eqs. (^Tj) and ( ^2| ) give different magnetoplasmon excitation branches: three 
triplet spin-density-excitations (denoted by ijJa± and lOz), and one singlet charge-density-excitation (denoted by ujp). 
In Fig. I we show the calculated dispersion energies of the charge mode, LOpiqi^), and the lowest energy triplet spin 
mode, LOaj_{q±), for a typical parallel magnetic field,i?|| — 11 Tesla at filling factor v = Q and other system parameters 
chosen to correspond to the experimental sample ||7||. The most important feature in the spectra is that there is an 
energy minimum (" magneto-roton" ) at a finite wavevector, q* ~ Z^^, in the spin mode dispersion along y direction 
(perpendicular to the in-plane magnetic field which is along the x axis), while no such finite wavevector minimum 
exists along the x direction. Comparing with the zero width 2D results (without any in-plane field) obtained in pO[ , 
where a roton- minimum is found in the spin mode along both directions, one finds that the finite width reduces the 
electron-hole binding energy (Eq. (|36|)), which is the origin of the roton-minimum in the magnetic exciton picture, 
along the direction of the in-plane magnetic field. (Note that for a ZWW system, the in-plane magnetic field does not 
change the electron orbital wavefunctions and it simply increases the Zeeman energy only, which is proportional to the 
total magnetic field) . From the energetic point of view, therefore, this "softening" associated with the development 
of the roton minimum (transverse to the in-plane field direction) implies that the ground state of such a quantum Hall 
system has a tendency to to make a transition from a uniform, unpolarized state to a spin density wave state with 
broken translational and spin symmetries, particularly if this roton-minimum reaches zero energy in some situations. 
In our calculation, the minimum energy of the spin mode (wo-^) goes to zero energy at B\\ — 12.5 Tesla. However, 
this value is close to, but slightly larger than the critical in-plane magnetic field, Bj| = 11.1 meV, where the ground 
state makes the first order spin polarization transition from a paramagnetic (z^| = v^) to the spin-polarized state 
(i/| = j/f — 2) (see Table ^ in the Hartree-Fock approximation. Therefore, within our HF approximation the roton- 
minimum of the spin mode dispersion does not actually go to zero energy before the whole system undergoes a first 
order phase transition to a polarized ground state. Calculating the collective mode energies for a polarized ground 
state after level crossing, we find that this roton minimum energy does not vanish, and in fact, may even increase 
in magnitude. Therefore we do not observe a true mode softening in the spin density excitation in the present 
Hartree-Fock approximation although we see a clear tendency toward such a possibility within our HF theory. It 
is certainly possible that a more sophisticated approximation going beyond the HF approximation would produce 
such mode softening (see the discussion in the following sections). Note that the charge collective mode energy in 
the long wavelength limit is exactly the same as the noninteracting energy separation, uj2 = 1-72 meV, in Fig. IJ, 
for results calculated in both the first order HF approximation and the SCHFA, within a 3% numerical error. As 
should be obvious from our results, there is no qualitative difference whatsoever between the results in these two 
approximations, which is not unexpected. Therefore, from now on, we will only show results obtained in the first 
order HF approximation, not only because of its compu tation al simplicity (saving considerable time in numerical 
calculations), but also because, as mentioned in Section III E| , we believe that the leading order HF calculation is 
really more consistent with our TDHFA theory for the collective modes. 

In Figs. H and ^ we show respectively the charge and spin mode dispersions for = 2, 4, 6, and 8 system by 
changing B± (total electron density is fixed) with all other parameters the same as in Fig. ^. The RPA peak is 
relatively weaker in stronger perpendicular magnetic field (smaller v) , while it is more pronounced when more Landau 
levels are occupied (larger v). On the other hand, the energy difference between the long wavelength limit (which is 
just the noninteracting energy gap, 102, according to Eq. ([5^)) and the roton minimum of the charge mode excitation is 
larger for smaller v (stronger B±) system. This indicates that the multiple absorption peaks observed in the polarized 
inelastic light scattering experiment should be separated more widely for smaller v (stronger B±). For the spin 
mode excitations shown in Fig. ^ the results for different filling factors are quite similar, except for their different q* 
(i.e. the position of the magneto-roton minima) due to different perpendicular magnetic field values. 

In Figs. and |^ we show respectively the charge and spin mode dispersions for ly — 6 system but with different 
confinement energy, cjq, as indicated in the figures. Larger confinement energies indicate smaller well widths in z 
direction. Therefore we have a continuous "transition" from 3D to 2D by increasing luq at a fixed density. This 
transition is observed clearly in Figs. |^ and ^ where the spectra in x and y directions become very similar for higher 
values of loq, reproducing the zero width (strictly 2D) results | [lO| . On the other hand, the roton minimum energy of 
the mode decreases for weaker confinement potential (larger effective well width), showing more of a tendency 
to have a spin-density-wave instability in a wider well. Another important feature can be seen in the charge mode 
dispersion. When the confinement potential is weak (e.g. luq — 7 nieV), the energy of the roton-minimum is smaller 
than the mode energy in the long wavelength limit {q± = 0). But the roton energy becomes larger than the long 
wavelength mode energy when the confinement potential is increased to = 19 meV, reproducing the results of the 
pure 2D system |10 , where the roton minimum is typically at a higher energy than the long wavelength mode energy. 
Therefore the finite width effect also enhances the tendency of a charge density wave instability against the ground 
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state. 

In Fig. ^ we show a typical singlet charge density magnetoplasmon mode (wp) dispersion of = 1 as an example 
of odd filling factors in TDHFA. For v = \, there will be no first order phase transition by Landau level crossing in 
any strength of in-plane magnetic field. When _B|| is more than 30 Tesla, we find a charge density wave instability at 
a finite wavevector perpendicular to the in-plane magnetic field. More detailed Hartree-Fock analysis shows that ]3^ ] 
this CDW is a kind of isospin skyrmion stripe, which has a charge density modulation in the x — y plane as discussed 
m Section 0. 

As a final remark, we note that Eq. (^9|) and Eq. (|4^) are based on TDHFA, which is exact only to the lowest 
order in the ratio of interaction energy to noninteracting energy gap {{e^ / cqIq) j 102) ■ Therefore it is a priori not clear if 
this leading-order many-body approximation can be used to study the mode softening phenomena near level crossing, 
where the interaction energy is necessarily comparable to (or stronger than) the noninteracting level separation since 
the noninteracting levels becomes degenerate at the critical point. However, to the best of our knowledge, no other 
systematic reliable technique is available to calculate the collective mode energy and such a mode softening behavior 
was earlier successfully treated within the TDHFA in the context of the second order phase transition related to 
the canted antiferromagnetic phase of a double layer system in the presence of interlayer tunneling and Zeeman 
splitting Q. Therefore we believe our results should be qualitatively valid in the level crossing regime. We do not, 
however, exclude the possibility that correctly including higher order interaction effects may very well reduce the 
roton-minimum energy zero at a finite wavevector before the system undergoes a first order phase transition to a 
polarized ground state. We do not know how to go beyond the TDHFA in a systematic current-conserving manner 
but speculate that such a calculation may very well give rise to a finite wavevector softening of the magneto-roton 
producing a quantum phase transition to the symmetry-broken phase. Our speculation is partly based on our finding 
that TDHFA actually predicts such a transition at B\\ = which happens to be sightly larger than the critical field 
for the first order transition. 



IV. SCREENING EFFECTS 



In the TDHFA shown in the previous sections, electron-electron interaction is the the bare Coulomb interaction 
without taking into account screening effects from the electron-hole fluctuations in the Landau levels. In this section, 
we will incorporate screening effects in our magnetoplasmon calculations. Actually, in Eq. (^9|) and (Q), a complete 
formula for the dielectric function in TDHFA has been given, but this formula is in general too complicated to be widely 
used in an integer quantum Hall system. In this section we will derive some convenient formulae for the dielectric 
function, e{q,uj), in different reasonable limits. Including such screening effects in the bare Coulomb interaction 
one may study the magnetoplasmon excitations beyond the time-dependent Hartree-Fock approximation, where the 
interaction used in the Green's function and vertex function is the unscreened one (see Figs. || and |o|(a),(b)). For 
convenience of discussion, we will first show the results for a zero width well, the system most theoretical researchers 
consider in the literature, and then the screening results for the WWW system of interest to us. It is shown below that 
for a ZWW, we can obtain a " scalar" (not matrix) dielectric function including both RPA and ladder diagrams shown 
in Fig. ^(c) in TDHFA, i.e. the screening effect in a ZWW is independent of the level index within TDHFA. This 
result is valid beyond the pure RPA result proposed before in Ref. [35| and should apply even at low density, where 
only a few Landau levels are occupied, because of the inclusion of the ladder diagrams (left out in Ref. ^^). For a 
WWW, instead of using the complete result shown in Eq. ( ^9| ) and (^^, we will derive a conventional formula in the 
strong parallel magnetic field region, which in some sense is effectively similar to a ZWW system as mentioned in Sec. 



II A , An analytical expression for the dielectric function can be obtained when only the RPA screening is considered 
(neglecting ladder diagrams) and is a good approximation for high density systems. Note that these general formula 
of screening effects could be used to study other interaction-induced electronic properties of quantum Hall systems 
1^ , ^ , and are therefore of broad general interest in quantum Hall problems transcending the specific applications 
we are dealing with in this paper. 



A. Screening in a Zero Width Well 



For a strictly 2D ZWW one can neglect the z degree of freedom completely, and therefore the interaction matrix 
element of Eq. (|2C|) and Eq. (pl|) can be simplified to the product of Coulomb interaction and the function (q) : 



U?i?ru:n.nJ.<l±) = {q^)Al^,J-q^)Al'^^M) = V'^{q)AZjq)Al'^^^{q), (57) 
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where V'^^{q±) = {2ttIqLz) ^ J dqzV{q±,qz) is the two dimensional Coulomb interaction and ^^{^ .(^l) = ^^f^ (?) 
is obtained by using the standard Landau level 2D single particle wavefunction in an integral similar to Eq. ([2]]). 

Its explicit formula can be obtained by taking the zero width limit [loq — > oo) of the function An^ujiq) in Appendix 
Note that in such a pure 2D system, electron wavefunctions obtained by SCHFA are exactly the same as the 
noninteracting wavefunctions, so that the results in the SCHFA and in the first order HFA are the same in this case. 
Wc use the superscript, "2D", to denote pure two-dimensional quantities in zero well width limit, and replace q± by 
its absolute value, q, in Eq. (|57|) due to the rotational symmetry in the x — y plane in the 2D limit. 
Applying Eq. (|57| ) to Eq. (^Sj) and Eq. ([39|), one obtains immediately 



U, 



2D 



{C0,q) ^UZ,^,,^nA^,q) ll-V'^{q) J] A^^^^Jq) ^ [f. 



r2D 



A. 



2D 



(58) 



where the matrix Y'^^{uj,q) is of the same form as Y(uj,q) in Eq. (^0[), but now with two-dimensional interaction 
matrix elements. The dielectric function for a ZWW system is therefore a scalar function and independent of the level 
index: 



e'''{0J,q)=ll-V'^iq) ^ A^^^^Jq) J] [f. 



^2D 



A 



2D 



(9) 



(59) 



and the corresponding irreducible polarizability, lif^{LO,q) can be easily obtained b y u sing e {uj,q) ~ 1 



V^^{q)nfX.{i^,q). Note that when the spin degree of freedom is considered in Eq. 



), all electron- hole pair 



fluctuations involved in e^^{uj,q) should be non-spin-flip pairs because Coulomb interaction does not flip electron 
spin. 

It is instructive to study the forms for the dielectric function in some special limits. First, in the low frequency 
region, where only fluctuations like {N, i) ^ {N + 1, |) and {N, 1) ^ {N + 1, t) are relevant, we can use the 2x2 
matrix of Fp(w, q) in Eq. (p2[) to express the dielectric function in the lowest order of (e^/eo^o)/w2: 



e^^iu; ^ 0,g) ^ 1 - V'^{q) \A%^^^,{q)\'Y.[Y^f{cu,q)]-' 



2V'^iq)\Aj,^^^,{q) 



AE 



ll,2D 
NN+1 



w^^(g) + iuj 
w^f (g) + icj' 



(60) 



where we have used the fact that Ai?^^'^^ = AE"^^'^^ for systems with even filling factors in the unpolarizcd ground 
state. 

Another good approximation for the dielectric function of Eq. ( ^9|) can be obtained in the high density limit, where 
it is well-known that the contribution of RPA diagrams dominates that of ladder diagrams in the correlation energy 
p7| . Starting from Eq. (|3^ ) and ( ^ and using iterations with Eq. ( ^ to represent „j_^„^ (g), which is now 

the same as J7^^„^,™„,„^(g), we have 
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(61) 



After retrieving the spin index, we obtain 
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= 1 + ^ (9)2.2. 2. {E,,„-E^„f + u^^\—) ' ^™ ' ^^'^ 

which is the same as the result in Ref. (using the identity: = (— l)'"(n!/(ri + m)!)a;'"L™(a;)), if we neglect 

the self-energy correction in the single particle energy, E^^^a (so that E^^^a = ^ = ("-c + l/2)w± — f'^z)- Note 
that the RPA result shown in Eq. ( |62|) includes the dressed single particle Green's function via the Fock self-energy 
correlation (the Hartree term is cancelled) , but it sums over all empty and filled levels and is therefore actually beyond 
the validity range of TDHFA which neglects multi-exciton effects. Both Eq. (|59| ) and Eq. (^) above are independent 
of the parallel (in-plane) magnetic field in the strict 2D limit, since the parallel magnetic field only affects the Zeeman 
energy in the strict 2D limit (and not any aspects of the orbital motion). However, the parallel field does, as expected, 
affect the dielectric function in a finite width well as shown below. 



B. Screening in a Wide Quantum Well 



For a WWW (specifically a parabolic WWW for our calculations), Eqs. ( p9[ ) and ( [40| ) show that the dielectric 
function is a matrix function strongly dependent on the level index of the interaction matrix element. In general, 
these expressions are not convenient for applications in different physical problems, and therefore we have to look 
for a good approximation for Eq. (|3^). First we could get a good low frequency approximation for the dielectric 
function by truncating the matrix size of Eq. (|39| ) into 2x2 and applying Yp of Eq. (^2|), i.e. considering electron- hole 



fluctuations only between the two nearest levels about the Fermi level. The result is similar to Eq. (60) 



eN.N+i{t^ ^ 0,qj_) ^ 7^-r——i 63) 



where LUa-Aqi_) and LUp{q±) are given by Eq. ( [15[ ) and Eq. (46) respectively. The difference between Eq. ( |60| ) for a 
ZWW system and Eq. (|6^ ) for a WWW system is that the former can be used for interaction between electrons in 
any Landau levels, while the later is correct only for electrons interacting between {N, "f (|)) and {N -I- 1, t (i)) levels 
in the low energy region of an unpolarized ground state. When considering higher energy excitation, say electrons 
from N — 2 level to level iV-|- 1, a larger matrix representation for the Yp matrix has to be used to get a self-consistent 
result, but it may exceed the validity region of TDHFA. It is instructive to check the asymptotic approximation of 
Eq. (|6^ ) in the static long wavelength limit by using Eqs. (p5|), (^6|), and (^): 

l^2 + ^N+l^^N^ -UN'N.N+lN+li^^^O) + '^UNN+l.N+lN{q±^0) 

eN,N+i['J,q± Uj ^ , , , ^HF w?F TTbmd 7T~^\ 

> 1. (64) 



UJ2 



^2 + ^N+1 - 



Note that eAr,Ar+i(0, ) does not go to unity because of the finite direct (Hartree) self-energy term, showing a 3D 
property. In a ZWW, however, the Hartree self-energy is a constant independent of the level index, and therefore is 
cancelled with each other in Eq. (|64[). When taking the large momentum limit {\q±\ oo), eArjv+i(0, gl) — *■ 1 for 

As in the ZWW, a scalar dielectric function similar to Eq. ( |59| ) can be obtained for a WWW system subject to 
a strong in-plane magnetic field. The similarity between these two systems is because the strong in-plane magnetic 



field effectively enhances the electron confinement energy of the well (note that ujb — y^o ~^ '^'\\ ^^"^ ^I^)- 
We start from the following general approximation (we use number labels (e.g. 1,2 • • •) and Greek labels (e.g. a, /3 
• • •) to replace the level indices, ni_2 -- and TOq,/?... for simplicity): 

U i4,Q/3 {q±)Up;,,23 iq± ) 

= J2^ Ai4(-'f±, -q'JV{q±,q'JAafi{q±,q'^) ^ Ap^{-q^, -q^_'')V{qj_, q^")A23{q±, q^") 



= J^YAi4{-q±,-qJV{q±,qJA23{q±,q'jYA^^{-qj_,-q^")V{q±,qz'')Aai3 x 



^23(gj-,gz") Aaf}{q±,q'^) 

z, . ^23(g±,90 ^a/3(g±,gz ) 

« ' 
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where C^^((fx, g^i Qz") has been approximated by its zeroth order value 



(65) 



(66) 



according to the exphcit expression of Aaf3 (q) shown in Appendix |^. Therefore the TDHFA screening similar to Eq. 
(p9|) for a WWW system could be obtained approximately as: 



e(cj,gl)- I 1 - XX [Fa/3,^i.(t^,gj_)] ^U^^^apiqi 

af3 fif 



(67) 



Similarly the high density approximation with RPA diagrams only can also be obtained by using the same approxi- 
mation: 
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(68) 



Comparing results of Eq. ( |62|) for a ZWW and Eq. ( |6q ) for a wide (parabolic) well, we find that the finite width effect 
enhances the anisotropy of the dielectric function through the coupling of x and z components of wavevectors. Note 
that Eqs. (|63|), (|67|), and (|6^) show no screening in the z direction because we have integrated out the z component 
in the interaction matrix clement by the single particle wavefunctions in Eq. ( |20| ) and have assumed the level index 
dependence of the dielectric function to be unimportant (see Eq. 66)). We believe that this is a good approximation 
for strong in-plane magnetic fields (see Fig. |l|), so that there is no appreciable static or dynamical polarization in the 
z direction to screen the Coulomb interaction. This approximation certainly fails when one wants to study excitations 
between levels of two different subbands in a weak in-plane field region. 



C. Numerical Results 



In this section, we show some numerical results of the collective mode energies including the screening effect. For 
convenience, we choose the dielectric function shown in Eq. (^) and consider static screening (w = 0) only. Therefore 
the algebraic matrix equations of Eqs. (|38|)-(^) are all of the same form except the Coulomb interaction is replaced 
by the screened one, V{q)/eRPA{q±,0)- However, the interaction of the RPA energy in Eq. (13^) and the Hartree 
self-energy are not screened in order to avoid double counting of bubble diagrams (see Figs. ]lO|(a) and (b)). We 
note that such screened TDHFA is not a strictly current-conserving approximation, because some other diagrams (for 
example, see Fig. |l^(c)) are not included, which may contribute to the same higher order effects as the screening 
bubbles. Therefore we can only estimate the screening effect to the magnetoplasmon energy qualitatively rather than 
quantitatively in our present study ||37[| . 

In the presence of screening, the first order phase transition point, B*^, moves higher values (see the fourth column 
of Table |) because the exchange interaction strength is reduced. This allows us to investigate the magnetoplasmon 
mode dispersion at higher values of in-plane magnetic field without changing the ground state configuration (i.e. 
avoiding the trivial first-order transition). In Fig. ^ we show the static dielectric function, e{uj — 0,q±) obtained 
by Eq. (|6^) in RPA for two different values of in-plane magnetic field at ^ = 6. For a stronger in-plane field, the 
screening effect is also stronger and more anisotropic. The anisotropic dielectric function shows that interaction along 
X direction (parallel to the in-plane field) is screened more than the interaction along y direction (perpendicular to 
the in-plane field). 



In Fig. 12 we show calculational results of the charge (ujp) and spin {uja+) mode magnetoplasmon dispersions 
including RPA screening effects (dashed lines) with other parameters the same as Fig. H For comparison, the 
unscreened results (dotted lines) and the screened results with higher in-plane magnetic field (solid lines) are shown 
together in the same figure. Comparing unscreened and screened results at =11 Tesla (dotted and dashed lines 
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respectively), one can find that the screening effect does lower magnetoplasmon energies for both charge and spin 
modes in the large \q±\ region due to the shrinking of Fock self-energy. But this effect is relatively weaker in the long 
wavelength limit (small \q±\ limit) due to the cancelation between the Fock self-energy and the electron- hole binding 
energy (the generalized Kohn's theorem). In the intermediate \q± \ region, the roton- minimum becomes less prominent 
than the unscreened result, and the dispersion becomes flat. Therefore, fixing all the other system parameters, the 
screening effect is not very important in determining the roton-minimum energy. On the other hand, as mentioned 
above, the screening effect reduces the electron self-energy and increases the critical value of for the unpolarized-to- 
polarized first order phase transition. Therefore one can, in the presence of screening, calculate the screened collective 
mode at higher in-plane magnetic field based on the same unpolarized ground state since the first order transition is 
now pushed to higher fields. In Fig. ^we show the result of magnetoplasmon dispersion calculated at = 12 Tesla 
(solid lines). The roton minimum of the spin collective mode {(J+) becomes lower than 0.1 meV, showing an almost 
mode softening at finite wavevector along the direction perpendicular to the in-plane magnetic field. Our results 
therefore indicate that inclusion of screening, effect, as well as lowering the confinement potential and increasing the 
electron density, could help to stabilize a new anisotropic ground state with broken translational and spin symmetries 
associated with the softening of the spin collective mode. Such a symmetry broken phase may very well be the cause 
for transport anisotropy observed in Ref. m . 



V. DISCUSSION 



In this section, we briefly discuss the possible phases of this new ground state based on our collective mode calculation 
results shown above. More detailed theoretical results on these exotic quantum phases will be given elsewhere [ p4[ . 
Similar to a DQW system ||2^], where the layer index is treated as an isospin degree of freedom, the level index of the 
closest two levels around the Fermi level can be used to construct an isospin (here it is also one-to-one related to a 
real spin) space, and create a coherent wavefunction for the possible new ground state in a single Slater determinant, 

I*,) = n( e^'^^^i COS ^4,,^^/,^ + e-^'^^^o sin ^4^,,^^^/,,, ) |0), (69) 

k 

where cj^ creates an electron in state 4>n,k,s{r) with spin s, and |0) denotes the ground state with + 1 filled 
Landau levels of spin up and N levels of spin down. We consider six different phases constructed from Eq. (|6S|), 
corresponding to different variational parameters, Wk and Q. When Wk is constant, the wavefunction of Eq. ( |69| ) can 
describe three non-stripe phases: (i) a fully (un)polarized uniform quantum Hall phases for Wk = (0)7r, (ii) a simple 
inter-level coherent phase for Q = and Wk 7^ 0, tt, and (iii) a spiral phase for finite Q and Wk 7^ 0, tt. When Wk 
changes periodically with fc, three different kinds of stripe phases arise: (i) simple stripe phase for Q = 0, which has 
no spiral structure, (ii) skyrmion stripe phase for finite Q, but Q J- h, where h is the normal vector of the stripe 
formation. Such a skyrmion stripe phase has both charge and spin modulation in different directions, and therefore 
has finite topological charge density oscillation in real space 1 26 3^,Q ; (iii) spiral stripe phase for finite Q with Q 



n. 



This spiral stripe phase has charge and spin modulation in the same direction but no topological charge oscillation 



in real space. We should point out that the wavefunction of Eq. (69) is based on a special choice of Landau gauge, 
A = {0,B±x — _B||Z,0), and therefore gives the stripe direction along y, i.e. perpendicular to the in-plane magnetic 
field. Choosing another kind of Landau gauge, where electron momentum is conserved along x, A ^ {—B±y, — _B||Z, 0), 
we can construct a stripe along x direction and the noninteracting Hamiltonian can be solved exactly by a canonical 
transformation ||3^. Then one can write the trial wavefunctions of these different phases and obtain their energies 
in Hartree-Fock approximation. The one of the lowest energy states should be the ground state near the degeneracy 
point, _B|| = _Bj|. Details will be presented elsewhere [Q. 

On the other hand, the magnetoplasmon excitation spectra we obtain in previous sections also gives us important 
information about the new ground state near the degeneracy point. First, the asymmetry of spin density mode 
in X and y direction and the near mode softening in y direction (shown in Fig. ^) strongly indicate that the new 
symmetry-broken ground state, if it exists, should have a spin spiral structure at finite wavevector in y direction. 
This may be a spiral spin density wave, when only one of the ordering wavevectors ±(0, q*) is present, or a coUinear 
spin density wave, when there is ordering at both wavevectors with equal amplitudes. The former can be visualized 
as a spin density wave where electron spin has a spiral structure around the total magnetic field direction in order to 
optimize the exchange energy. Therefore collinear spin density wave, spiral, skyrmion stripe, and spiral stripe phases 
are the possible candidates for the symmetry-broken phase. As for the existence of any possible charge density wave 
instability, we could not obtain much information from our collective mode calculation in TDHFA. But it is apparently 
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true that interaction effects are more important for stronger in-plane magnetic fields, where the noninteracting energy 
separation, W2, becomes very small. 

Considering the experimental results Q , where the resistance along the in-plane magnetic field becomes finite when 
the in-plane magnetic field exceeds a critical value, we find that the stripe formation, if it exists, should be along the 
direction perpendicular to the in-plane magnetic field (i.e. its normal wavevector is in x direction) to produce such a 
transport anisotropy. Using the above result that the spin density modulation has a wavevector, q*, in y direction, we 
find that only the skyrmion stripe phase is consistent with all of these constraints and should be the best candidate 
for the new ground state. Although our Hartree-Fock calculation shows that the spiral phase has slightly lower energy 
than the skyrmion stripe phase ||2^ , we believe that this may be due to the nonparabolicity of a realistic WWW or 
the correlation effects not included in our HF approximation. We therefore speculate that the anisotropic ground 
state observed in Ref. is our proposed novel spin skyrmion stripe phase. This may also be true for the transport 
anisotropy earlier observed in Si based 2D systems, but the additional complications of valley degeneracy in Si 
makes the application of our theory mode different. 



VI. SUMMARY 



We study the magnetoplasmon excitations of a parabolic quantum well system in a tilted magnetic field. Starting 
from the many-body theory in coordinate space, we integrate out the continuous variable and obtain an algebraic 
matrix representation of the dielectric function and hence the magnetoplasmon mode dispersion in TDHFA. Focusing 
on even filling factors, a roton-minimum near zero energy in the spin channel is observed at finite wavevector along the 
direction perpendicular to the in-plane magnetic field. By changing the confinement potential, we have a continuous 
transition from a 3D plasmon excitation to the pure 2D results in our calculation. Including the screening effect, 
which is another important part of our work, we find that the roton-minimum energy could be even more suppressed. 
Although it does not reach zero energy before possibly undergoing a first order phase transition from an unpolarized 
ground state to a polarized one, its small excitation energy at finite wavevector suggests a possible spin-density- 
instability to an exotic symmetry-broken ground state in realistic systems. We discuss various phases that may result 
and propose that the recent transport anisotropy measurement in experiments [Q can be explained by a skyrmion 
stripe phase, where spin and charge density modulations are in different directions. The theoretical technique used 
in this paper could also be used to study other quantum Hall systems in quasi-2D quantum well nanostructures. In 
particular, our screening theory is more complete than the existing theory, and should have wide applicability. Finally 
we point out that our predicted collective mode dispersion may be directly verified via the inelastic light scattering 
spectroscopy. 
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APPENDIX A: SELF-CONSISTENT HARTREE-FOCK EQUATIONS 

In this section we derive the self-consistent Hartree-Fock equations for the single particle wavefunctions. Starting 
fromEq. @, we can use a Fourier transform of V{f) to obtain 



En,a4>n,k,cr{r) 



q fh.p.s 
q m,p 



<i}n,k,a{r) 



4'n,k,cr{r) 
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(Al) 



where the form function, A'^^{q), has been defined in Eq. (pl|). Assuming a uniform ground state, we can separate 

4>n,k,iy{'r) into a product of a plane wave, e*'^^/-\/X^, and the hmction, '^fi^„{x + klQ,z), which satisfies the following 
eigenvalue equation: 



no + E ^fe) E '^^x^'A<i{q.) e'^^' 



where Vf^^^ is the filling factor of Landau level fh and spin cr, satisfying 



(A2) 



(A3) 



to conserve the total electron density. 

Now we expand z) in terms of noninteracting wavefunctions with the same spin (note that V(r ) allows no 

spin- flip, so that a is conserved and no spin hybridization occurs): 



where | • • •)o represents a noninteracting eigenstate. We have 

= E E^"^' ^1^1' o(^2, <j\n, cr)yl?^:J'^(q ). 
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Using Eq. ( [Aq ) and multiplying by the noninteracting wavefunction from the left of Eq. (A2), we have the self- 
consistent Hartree-Fock equation in a matrix representation: 
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(A6) 



where we have used Eq. ( pO| ) and Eq. ( ^l| ) to express the direct and the exchange potential. Eq. ( [A^ ) is the matrix 
representation of the Hartree-Fock Hamiltonian in our system, which should be solved self-consistently to get the 

energy eigenstate via vector elements, o{rn, a\n, a). 

Another expression for the eigenenergies can be obtained directly from Eq. (A2), by integrating another eigenket 
^n,a{x, z) from the left. We obtain 



irT.^(i^) E ^™,.'A^:;^:(g.)A^~^(-g.)- ^E^(9)E^™^-^;?rn(M^^ 



(A7) 



21 



where UZ%.^ .{6^) and u'J^^Z'^^iO ^) are those defined in Eq. and Eq. (^. 

Using this self-consistent Hartree-Fock equation, Eq. (A6), it is also easy to include any nonparabolic effects of the 
realistic confinement potential, U{z). Assuming the deviation of the realistic U{z) from a parabolic one, Up{z), to be 
small, i.e. \AU{z) — U{z) — Up{z)\ <C loq, we can calculate its matrix element. 



{n' ,a\/^U{z)\m' ,a) = j dx j d^; $s/ (x, z) AC/(z)$^ / (x, z) 

= j dx j dz^n,{x,z) -^^Af/(q,)e^' 



^A'{X,Z) 



^ ^ ^ 1 /" 

= —Y,^U{q,)Aft'^'iO±,-q,)^— / dq,AU{q^)Afi'^'iO^,-q,). 



(A8) 



and incorporate it in Eq. (A6) to calculate the self-consistent Hartree-Fock eigenenergies and eigenfunctions. In all 
our numerical work presented in this paper, however, we have taken U{z) to be parabolic throughout. 



APPENDIX B: MAGNETOPLASMON EXCITATION ENERGY THROUGH THE MAGNETIC 

EXCITON WAVEFUNCTION 

In this section we show that the magnetoplasmon excitation energies both in a thin 2D (ZWW) well in only a 
perpendicular magnetic field (situation discussed in |[lO|) and in a wide parabolic well with a tilted magnetic field 
(situation discussed in this paper) can be written in a simple and instructive form by using exciton wavefunctions 
proposed in and its appropriate WWW generalization constructed in our Eq. (p^), respectivel y. For the first 
case, we take the static exciton wavefunction suggested by Kallin and Halpcrin in Eq. (2.9) of Ref. and set the 
center of mass coordinate and the total momentum of excitons to be zero: 



vI/2^,„JAx,Ay)^ / dr;e-''^^/'o^(°)(ry-f Aa:/2X)(,y-Aa:/2), 



(Bl) 



where Aa; and Ay are the relative coordinates between the hole in a filled level (denoted by — n) and the electron 

in an empty level (denoted hy np — na + m) ; (x) is the wavefunction of one-dimensional single harmonic oscillator 
as shown in Eq. (@) with li replaced by Iq. In the lowest order of (e^ /dQ)/u)±, there are four distinct contribution to 
the magnetoplasmon excitation energies: noninteracting energy separation, exciton binding energy, RPA energy, and 
exchange self-energy [ p^ , 



,2D 



^E^Uq) 



A rpmn 



(B2) 



where the last three terms can be re-expressed in terms of ^' 



2D 



as follows 



dArj_V'^"{Ar- l'^qi_ x z) \^ 



2V^"{q) 
27r?2 



2D 

n+m.n 



exch n+m 



- = 



27r/2 



dArj_V'iAf±_) 



^2D 

n-\-m,n-\-m 



-<„(Af_,) E *?.f*(Af^) 



(B3) 
(B4) 



(Afj_) 



(B5) 



where N„^^^^ is the level index of the highest occupied Landau level with spin aa{fi)- Interpretation of the formulas 
in ( p33| ) and ( B4 ) is straightforward. The binding energy integrates over relative positions of electron and hole in the 
exciton, whereas the RPA term involves electron and hole annihilating each other and is proportional to the probability 
of finding two particles at the same position. Ai?™"^ in Eq. (85) is the difference of exchange self-energies between 



the two relevant levels, and indicates the relative many-body level shift. The exchange self-energy of level n, 
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expressed in Eq. ( B5 ) can be understood as the integral over relative positions of electrons between level n and the 
lower levels, I, of the same spin. Note that uj'^{q± = 0) = muj± for m = 1 in the charge mode channel, satisfying 
Kohn's theorem for this ZWW system. The equivalence between above expressions and the results in can be 
easily seen by direct substitution. 

For a parabolic well, the magnctoplasmon energy expressed by the magnetic exciton wavefunction (sec Eq. ([l^)) 
can be obtained by using similar notations as above (let ria — n and np = n + m to denote the hole and electron level 
indices): 



mn 
direct 



where AE^ 



iiLOi + TO2^2 — ^zicrp — <7a) is the uoninteracting energy gap between the two levels, and 



AEZSi'l^) 



AE-liq^) 



21111 



dArViAr-llq^x z) / dZ l^ftn.nJAx, Ay, Z, Az)\ 



dAfT/(Af) e'('-^-+«='^^) j dZ^SJn,,nJ-qyllqJlZ + Az/2,0) 



>^n,,nA-Qvll<lJlZ - Az/2,0) 
-1 



(B6) 
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\ -Tirana 
direct 



-n^^n^ {Ax, Ay, Z, Az) ^ vj/-^ (Ax, Ay, Z, Az) 

la 

^ I dAfF(Ax, Ay, Az) V ( dZ*;-;-(0, 0, Z - Az/2, 0) 

T 

X (0, 0,Z + Az/2, 0) - (0, 0,Z + Az/2, 0)] , 



(B9) 



(BIO) 



where the summation over I means the summation over all occupied levels with quantum number, (^1,^2), and sum- 
mation over la{p) is the summation of all occupied levels with the same spin as the state ria^py The interpretation of 
these equations is simi lar to th e ze ro wid th si tuation, except for an extra integration over z coordinates. 

Note that EqS;_^(|B|), (||), (^) and ( plO| ) can be transformed to the momentum space by using the A function 
defined in Eq. (|l|): 



Anpiia {q) 



dx J dze-''^-^~''^^'^ft,{x - qyPo/2,z)^ftAx + qyll/2,z) 
dz e""^"^5's^, {-ly^l, Qxll.z, 0), 



(Bll) 



so that we obtain 
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2ttI^L^ 



^V{(lx,qy,Pz)\ Afip n^{qx,qy, Pz 



ae: 



direct 
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0^2 



which are identical to the results we have derived before in section [[III] by noting that At^^^q) — A^^fi {—q). 
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APPENDIX C: ANALYTICAL EXPRESSION FOR A^°'>'''S' (Q) 

IV_ w„ ' 



The explicit formula for the function J^^'Z" [q) we use in this paper can be evaluated by using the known math- 
ematical properties of the generalized Laguerre polynomial. Since it is defined by the noninteracting wavefunctions, 
which are not dependent on the spin index explicitly, we can neglect the spin index totally here and calculate the 
orbital integer directly from Eq. (|47|). Using Eq. and Eq. (^), we obtain the following results (for convenience, 
let fia = {na,n'^) and — {np.n'p)): 



I 



X exp 



4Ai 



exp 



sin 6{qylaY + {smOq^lo + cosOq^h) \2 
4A^ 



=F coa9{qylo) — i{cos9qxlo — sin9qzlo)Xi\"'""^ f =F sin6{qylo) — i{sm9qxlo + cos 6'gz^o)A2 V""" 



/2AT 



2 nf„ 1 ^2 



cos^6'(%/o) + icos9q.Jo ~ sin 9q Jo) Xi\ j^rn'^^ ( sin 9{qylo) + {sm9qJo + cos9qJoyX 



, (CI) 



2Ai J V 2A2 

where ± is the sign of n^'' — ni'^ for each bracket and n^'l, ■ , n = Min{Max){n}^ ^ n'R^-, and m*-'), = \n'a — rb-a\- 

-^1,2 = ('1,2/^0)^ are dimensionless parameters. L^(x) is the generalized Laguerre polynomial. 
As for a ZWW, we can let ujq ^ 00 and obtain 



exp 



„2 72 1 

4 



±qyla - iqJo 
V2 



1 '0 

2 



(C2) 



where q= \q±\^ and all notations are the same as in Eq. (CI) above 
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v 


noninteracting 


interacting 


interacting 






(unscreened) 


(screened) 


6 


19.8 


11.1 


12.2 


8 


19.8 


10.4 


11.5 



TABLE I. Table of the critical values of the parallel magnetic field, Sj| , where a first order phase transition occurs from an 
unpolarized ground state to a polarized one for the parameters of Ref. [7] . 
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FIG. 1. Calculated landau level energy spectra for noninteracting electrons in a parabolic quantum well with a parallel 
(in-plane) magnetic field, B\\ . The system parameters are chosen to be the same as the experimental data in ref. [7] for = 6. 




FIG. 2. Feynman diagrams of the time-dependent Hartree-Fock approximation. Solid lines are single particle Green's 
function and wavy lines are Coulomb interaction. Single(double) lines are bare(dressed) Green's function and/or interaction: 
(a) the self-consistent Hartree-Fock approximation for the single electron Green's function; (b) and (c) are respectively the 
Dyson's equations for electron-electron interaction, single electron Green's function and vertex function in the time-dependent 
Hartree-Fock approximation. The second term of (c) is the ladder series, while the third term is the bubble series (RPA diagram), 
which does not appear when calculating the vertex function for spin-flip excitations (since the interaction is spin-conserving) 
as mentioned in the text, (d) is the Green's function in the first order Hartree-Fock approximation. 
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FIG. 3. Energy level configuration for electron-hole pair excitations. Solid(dashed) lines are for spin down(up) levels with 
level index in the left hand side (the first orbital level index is set to be zero) , and the upward arrow represents an electron- hole 
excitation (a magnetic exciton). (a) is for the two 2x2 matrix representation of Eqs. (^l|) and ji^): electron-hole pairs of 
numbers 1 and 2 are for Ya, and numbers 3 and 4 are for Yp respectively, (b) shows the configuration for one spin-flip excitation 
{Sa = +1) including next higher order energy excitations, which are beyond the TDHFA developed in the paper. (Note that 
in (b), the excitation from level n to n + 2 does not couple to pair number 1 due to parity symmetry in a parabolic well.) 
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FIG. 4. Magnetoplasmon dispersions for u = 6 and By (along x direction) is 11 Tesla, calculated from Eqs. (^) and (^6|). 
(a) and (b) are for momentum along x and y directions respectively. Thick(thin) lines are for wavefunctions calculated from 
self-consistent Hartree-Fock and from first order Hartree-Fock approximations respectively. 



29 




( 1 cm 



; '(b)' q, = C 










i/ = 6 1 















(10** cm"') 

FIG. 5. Charge mode dispersion, tjp((fx), of magnetoplasmon excitations of the same system as used in Fig. ^ but for 
different filling factors, v = 2, A, & and 8, for comparison. 
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FIG. 6. Same as Fig. ^, but for spin mode, u}a^{q±_), dispersion. 
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FIG. 7. Charge mode dispersion of magnetoplasmon excitation of the same system as used in Fig. W but with different 
confinement energy, loo, at filling factor, ly = 6. Zero field well widths are about 260, 200, 175, and 155A, corresponding to 
LJo — 7, 11, 15, and 19 meV respectively. The parallel magnetic field is 11 Tesla for all results. 
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FIG. 8. Same as Fig. |^, but for spin mode dispersion. 
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FIG. 9. Charge mode magnetoplasmon dispersions for u — 1, B_l=3 Tesla, and tJo = 3 meV, calculated in the TDHFA. 
Solid, dotted, and dashed lines are for _B|| = 20, 25, and 30 Tesla respectively, showing a charge mode softening in y direction, 
perpendicular to the By direction, (a) and (b) are for wavevectors along x and y directions respectively. 




FIG. 10. (a) The screened exchanged energy and (b) the screened ladder diagrams used in the screened TDHFA developed 
in Sec. IV (see Fig. The interaction lines of the direct (Hartree) energy and the RPA diagrams are not screened to avoid 



double counting, (c) A diagram not included in the screened TDHFA but of the same order as a screened ladder diagram shown 
in (b). All notations are the same as those in Fig. 
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q ( ' 0*" cm"' ) 

FIG. 11. Static dielectric function in momentum space for v = 6 and cjo ~ 7 meV. Solid and dashed lines represent 
^{<lx,qy = 0) and e{qx = O,^^) respectively. Thick and thin lines are for = 11 and 8 T. 




FIG. 12. Dispersions of magnetoplasmon excitations for = 6 in both charge (p) and spin ((J+) modes including RPA 
screening (Eq. (pSl)) of the Coulomb interaction for _B|| = 11 and 12 Tesla (dashed and solid lines respectively). Results of 
unscreened dispersion are also shown (dotted lines, the same as Fig. ^)) for comparison. 
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